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Abstract 

We derive the scattering amplitude for Goldstone bosons of chiral symmetry off the 
pseudoscalar charmed mesons up to leading one-loop order in a covariant chiral effective 
field theory, using the so-called extended-on-mass-shell renormalization scheme. Then 
we use unitarized chiral perturbation theory to fit to the available lattice data of the 
S-wa.ve scattering lengths. The lattice data are well described. However, most of the low- 
energy constants determined from the fit bear large uncertainties. Lattice simulations in 
more channels are necessary to pin down these values which can then be used to make 
predictions in other processes related by chiral and heavy quark symmetries. 
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1 Introduction 


The hadronic interaction between charmed D mesons and the Goldstone bosons (p of the 
spontaneous breaking of chiral symmetry of the strong interaction interaction for short 

hereafter) is important for the understanding of the chiral dynamics of quantum chromody¬ 
namics (QCD) and the interpretation of the hadron spectrum in the heavy hadron sector. 
Many investigations have been devoted to study it in the last decade, partly triggered by 
the observation of the charm-strange meson T)*q( 2317) with = 0+ in 2003 [Il[2]. The 
T1 *q( 2317) couples to the DK channel, and being below the DK threshold it decays into 
the isospin breaking channel DsTt. In order to unravel its nature, theorists study the D-cp 
interaction and intend to extract the information encoded in it. For instance, the ZI*q( 2317) 
is interpreted as a DK molecule [3] by using a chiral unitary approach to the S-wave D-(p 
interaction In these works, the leading order (LO) amplitudes from the heavy meson 

chiral perturbation theory (ChPT) [THS] are used as the kernels of resummed amplitudes. 
Extensions to the next-to-leading order (NLO) can be found in Refs. |initl4| . 

Recently, renewed interest was stimulated due to the occurrence of lattice QCD calcula¬ 
tions of the scattering lengths given in Refs. |15)I16|. In these works, only channels free of 
disconnected Wick contractions are calculated, which are Dtt with isospin I = 3/2, DK with 
1 = 0 and 1, DgK and DgTT. There have been lattice results on channels with disconnected 
Wick contractions, such as Dtt with 1 = 1/2 [T7] and DK with I = om- With these lattice 
calculations, more insights were gained into the nature of the Z)*q(2317). The DK isoscalar 
scattering length was calculated indirectly in Ref. m, which is consistent with the result 
from the direct lattice calculation in Refs. [I8l[l9]. A reanalysis of the lattice energy levels 
for the D^*^K lattice data m was performed in Ref. [20] in terms of an auxiliary potential 
and an extended Liischer formula. These results suggests that the Z)*q(2317) is dominantly a 
DK hadronic molecule. 0 

The lattice data can be used to determine the low-energy constants (LECs) in the chiral 
Lagrangian of higher orders. Especially, the lattice data in Refs. [IMI6| were used in Refs. \P2\ 
dSHMSTj. In the majority of those investigations, unitarized extensions of ChPT, see e.g. 
Refs. [28p29j . are adopted so that one can consider larger meson masses and channel couplings. 
The unitarized chiral perturbation theory (UChPT) is especially necessary for the chiral 
extrapolation of scattering lengths in question since for larger quark (or meson) masses the 
interaction normally becomes stronger and could even be nonperturbative. However, in all 
the calculations in the framework of unitarized ChPT, the kernel of the resummed amplitude 
was only calculated up to NLO at most and is purely tree-level. Here, we will extend the 
calculation to the leading one-loop order, which is the next-to-next-to-leading order (NNLO). 

It is well-known that ChPT [3nH32j has become an useful and standard tool in study¬ 
ing the hadron interaction at low energies. Based on Weinberg’s power counting rules m. 
great achievements have been obtained both in the pure mesonic sector and the one in¬ 
cluding matter fields such as baryons, the latter known as baryon ChPT. There is a no¬ 
table power counting breaking (PCB) issue in baryon ChPT [33]: using the dimensional 
regularization with the modified minimal subtraction (MS) scheme in calculating loop inte¬ 
grals, the naive power counting does not work and all loop diagrams start contributing at 
O (p^), with p being a small momentum. There have been several solutions to this problem; 

^ A method of extracting the probability of a physical state to be a hadronic molecule in lattice using twisted 
boundary conditions is discussed in Ref. m where the H*o(2317) is used as an example for illustration. The 
H*o(2317) was also studied in a finite volume in Ref. |22| . 
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heavy baryon (HB) approach [341135] . infrared regularizaion (IR) |36j and extended-on-mass- 
shell (EOMS) scheme |37j (for a review and a detailed comparison of these approaches, see 
Ref. [55]h 

Likewise, ChPT including the heavy D mesons encounters the same PCB problem. To 
remedy it, in Ref. [23], the D-(j) scattering lengths were calculated in the framework of non- 
relativistic heavy meson ChPT [7H9| to the leading one-loop order in the heavy quark limit. 
Nevertheless, as mentioned by the authors and confirmed by Ref. [24], this nonrelavistic for¬ 
mulation neglects sizable recoil correction^. The calculations in various unitarized versions 
of ChPT in Refs. |121ll6p251l26j are performed in a covariant formalism, but only up to NLO as 
mentioned above. The first NNLO calculation of the scattering lengths was given by Ref. [23] 
using the EOMS scheme. However, the calculation in that work is perturbative while the 
interactions in certain channels are definitely nonperturbative. For instance, in the channel 
with {S,I) = (1,0), where S and I represent strangeness and isospin, respectively, the ex¬ 
istence of the Z)*q( 2317) below the DK threshold calls for a nonperturbative treatment of 
the DK interaction or inclusion of an explicit field for the I1 *q( 2317). In addition, all the 
NNLO counterterm contributions are neglected in Ref. |23| due to the poorly known LECs. 
In this paper, we intend to present a detailed covariant description of the D-(j) interaction up 
to NNLO in the framework of UChPT, and the EOMS approach which preserves the proper 
analytic structure of the amplitudes will be used in renormalization procedure. 

First, we will calculate the D-cf) scattering amplitude in covariant ChPT up to the NNLO. 
To our knowledge, the D-(j) scattering amplitudes (without vector charmed mesons) shown 
in the present work are the first analytical and complete results up to NNLoH The vector 
charmed meson contributions, surviving in the heavy quark limit, are also taken into account 
numerically to estimate their influences, although it was shown in Ref. m that their con¬ 
tribution to the 5-wave scattering is small. Renormalization will be performed using the 
EOMS scheme and it will be shown explicitly that the UV divergences are cancelled properly 
and PCB terms are absorbed exactly by the counterterms, which ensures that the EOMS- 
renormalized D-cp scattering amplitudes possess the proper analytic, power counting, and 
scale-independent properties. 

We will fix the values of the LECs by fitting to the available lattice data of the 5-wave 
D-cj) scattering lengths. Since the lattice calculations are performed at large unphysical quark 
masses, the perturbative expansion to a certain order may fail to converge. One way to 
solve this issue is to employ unitarized amplitudes instead of the perturbative ones. Many 
unitarization methods have been proposed in the past. In Ref. [29|, a unitarization approach is 
developed and used to study the KN interaction. The unitarized amplitude can be matched to 
the perturbative amplitude order by order. Throughout this paper, we will call this approach 
UChPT for convenience. In Refs. msi, the inverse amplitude method (lAM) was proposed 
and adopted to study the Trvr and Kn scattering. For the purpose of comparison, both the 
two approaches will be employed. 

This paper is organized as follows. In Section O the power counting and and its breaking 
by heavy meson masses will be explained briefly and the chiral effective Lagrangian will be 
given up to NNLO. In Section [S] details on the computation of the D-cp scattering amplitude 

^It is, however, known since a long time that in the heavy baryon approach such recoil corrections can 
easily be incorporated by using as the propagator i/{v ■ I — r /2m) instead of simply i/v ■ I [39] 

®The analytical expressions for the amplitudes involving vector charmed mesons, which survive in the heavy 
quark limit, are too lengthy to be shown explicitly in the paper and can be made available upon request from 
the authors. 
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using the EOMS scheme are exhibited. Together with the loop results shown in Appendix iBl 
the minimal but complete set of scattering amplitudes are given explicitly. Section 0] discusses 
how to obtain the partial wave amplitudes with definite strangeness S and isospin I from the 
physical process amplitudes, and the two unitarization approaches mentioned above will be 
introduced. In Section 01 the S'-wave scattering lengths are calculated and htted to the 
available lattice data at a few values of the pion masses, and the contributions of the vector 
charmed mesons will also be discussed. Finally, Section [6] comprises a summary and outlook. 
Some technicalities are relegated to the appendices. 


2 Theoretical framework 


2.1 Power counting and power counting breaking terms 


We denote the D-cj) interaction as Di{pi)(j)i{p 2 ) —)• D2{p3)4>2{P4:)- The scattering process 
is on-shell, hence, pj = M|,^, p^ = M|^, pj = and with (M^i) and 

M/jj being the masses of the incoming and outgoing D mesons (Goldstone bosons), 

respectively. In addition, the Mandelstam variables are defined as 

•s = (Pi + P 2 f , t = {pi- , u = {pi- p^f , (1) 


which satisfy the relation s -\-1 + u = . At low energies, one has 
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where A^ 


X 

{47rF,r) ^ 011 ^ 02 } denotes the high energy scale, with the pion decay constant 
Fjr ~ 92.2 MeV. The above small quantities can be simultaneously adopted as expansion 
parameters. In a more conventional notation, one denotes the small parameters by a unique 
symbol, say p, so that the power counting rules for the basic quantities read 


Mdi ~ 0 (p°), Md2 ~ 0(P°), ~ 0{p^), ~ 0{p^), t ~ O(p^), 

s - ^ 0{p^), s - 0{p^), u-0{p^), u-~ 0(p^)- (3) 


It is worth noting that the the chiral limit masses of the charmed mesons are of the same 
order as the corresponding physical masses. Every physical observable therefore has its own 
chiral dimension by using the above given power counting rules. 

Furthermore, in ChPT a power counting rule is assigned for each Feynman graph. In the 
present case, the chiral dimension n for a given graph can be evaluated from 


n = iL + ^ 14 — 2/0 — Id , 
k 


(4) 


where L, I4, /0 and Id are the numbers of loops, order vertices, Goldstone boson propa¬ 
gators and charmed meson propagators, respectively. 

For a specific physical observable, if there exist terms with chiral dimensions obtained 
using Eq. ([3|) lower than that given by Eq. ([4]), those terms are called PGB terms. The PGB 
terms show up only when there are heavy particles with nonvanishing chiral limit masses in 
loops as internal propagators. In our present calculation, since the heavy charmed mesons 
are involved in some of one-loop graphs, there will be PGB terms if we use dimensional 
regularization with the MS scheme. These terms can be treated in the so-called EOMS 
scheme, which has a power counting consistent with Eq. Q, as will be detailed in Section 01 
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2.2 Chiral effective Lagrangian 

The pseudoscalar charmed mesons can be collected in a SU(3) triplet, D = {D^, D~^, D^), 
and the light Goldstone bosons are in an octet, 



(>“ + 75’' 

7r~^ 
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TT~ 

1 

O 

+ 



1 

>0 



The chiral effective Lagrangian for D-<j) scattering can be decomposed into D meson-Goldstone 
boson interacting parts and pure Goldstone bosonic parts, which has the following form: 

i-eS - + ^D4> + ^00 + ^00 + • • • • (Oj 

Here, the numbers in the superscripts stand for the chiral dimensions, and the ellipsis denotes 
the higher-order chiral operators which will not be used here. Besides, the operators with 
external fields are also dropped (except for the scalar external field which is used for the light 
quark mass insertions). 

The familiar lowest order chiral Lagrangian for the Goldstone boson sector reads 

41 = T ^ (xt/t + Ux^) , (7) 

with U = exp (z\/2</'/To) and x = 2i?o diag(mu, rris). Here (...) denotes the trace in 
the light-flavor space, Fq is the pion decay constant in the chiral limit, and Bq is a constant 
related to the quark condensate. We will work in the isospin limit with and neglect 

the electromagnetic contributions. 

The 0{p^) pure Goldstone boson Lagrangian is needed for renormalization. Its LEGs 
enter the D-cj) amplitudes merely through the wave renormalization constants and the decay 
constants of the Goldstone bosons, which can be found elsewhere, see e.g. Ref. [32]. The 
relevant terms read 

^00 = ^4 (xf/t + Ux^) + Ls (5^G(5^G)t (xt/^ + Ux^) ) + .... ( 8 ) 

For the interaction in the D meson-Goldstone boson sector, the LO effective Lagrangian 
takes the form 


= V^DV^^D^ - M^DD^ , (9) 

where Mq is the mass of the D mesons in the chiral limit, and the covariant derivative acting 
on the D mesons is defined by 

V^D = D{d^ + Tt) , = {d^ + T^)D^ , (10) 
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with the so-called chiral connection -|- udfj_u^) /2. The NLO Lagrangian reads [TT|R 

^ i-ho{x+) - hix+ + h2{Uf,u^) - hsUf,u^^) 

+V^D {hiiu^u’') - h5{u^^, u"}) , (11) 

where the building blocks of the chiral effective Lagrangian are given by 


Un = i 


i — u 


u = exp 


( 

\V2Fo 


X± = u^XU^ ± uxu 


( 12 ) 


Here, the definition for x- is also given as it is needed for the NNLO Lagrangian which, 
following the procedure detailed in Ref. [32], can be constructed as 


£(3) _ £) 

^D4> — ^ 




igi[X-,Uu]+ 92 {[u^, ['Fu,uf^]] + [Uf,, [Vf^,u^]]) 
+g3D[Uf,,[V,,Up]]V>^'^PD^ , (13) 

where the totally symmetrized product of three covariant derivatives is defined as = 

{Fp, {F>u, 'Dp}}. 


3 D-(f) scattering amplitudes up to NNLO 

In this section, we exhibit the complete set of independent D-4> scattering amplitudes on 
the basis of the physical states. They correspond to 10 physical processes as listed in the 
second column in Table [H All the other amplitudes can be obtained by using either crossing 
symmetry or time-reversal invariance. In what follows, we will hrst calculate the tree-level 
amplitude which can be reduced into a common structure but with different coefficients 
because of SU(3) flavor symmetry. Then the loop amplitudes will be given explicitly. In the 
end, the renormalization procedure within the EOMS scheme will be discussed. 

3.1 Tree-level contribution 

The Feynman diagrams of the tree-level contribution to the scattering amplitudes are 
displayed in the first line of Fig. 13.11 Since we do not consider the exchange of resonances, such 
contributions are encoded in the contact terms for the D-cj) scattering. When calculating the 
Feynman diagrams, all the bare parameters, such as the decay constant Fq and the masses, 
are maintained. They will be replaced by the corresponding physical quantities when the 
renormalization is performed. The LO, i.e. 0{p), tree amplitude is the Weinberg-Tomozawa 
termEI, and has the following form, 

A^^\s,t,u) =Clo^^ , (14) 

'^As in Ref. [1^, the he term in Ref. m is dropped, and the x+ = X+ ~ (x+)/3 is replaced by x+ which 
amounts to a redefinition of ho and hi. The he term is redundant, since 

— heT>fj.D[u'^ ,u'YDi,D^ = ^ ^D[u'^ -h -I-higher order terms, 

where the first term is zero due to the symmetry property of the Lorentz indices g , v, and the higher order 
terms are contained in the higher order Lagrangians. 

®As the vector charmed mesons are not taken into account, there is no Born term due to the exchange of 
these mesons. 
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Tree (1) 


Loop (1) 


Tree (2) Tree (3) 



Loop (2) Loop (3) 


Loop (4) 


Figure 1: The 1-point irreducible (IPI) Feynman diagrams for D-cj) scattering up to leading 

one-loop order. The solid (dashed) lines represent the D (Goldstone) mesons. The square 

( 2 ) 

stands for the contact vertex coming from Lagrangian while the filled circle denotes an 
insertion from All other vertices are generated either by or ^ 00 - 




where the coefficients Clo different physical processes are listed in Table [TJ The Weinberg- 
Tomozawa term depends only on the pion decay constant due to the fact that it originates 
from the kinetic term in 13^0, which is a result of the spontaneous breaking of chiral symmetry 
in QCD. 

The 0{p^) Lagrangian T ^0 generates the tree-level contribution at NLO as 


To L 


(15) 


where the coefficients are shown in Tabled] and the functions and are 

defined by 


H 2 i{s, t, u) = 2 /l 2 P 2 • P4 + hi {pi ■ P2P3 ■ Pi + Pl ■ PiP2 ' Ps) , 
H 35 {s, t, u) = h3P2-pi + hs (pi • P2P3 ' P4 + Pi ' PiP 2 ' Ps) ■ 


Finally, the tree-level amplitude at 0{p^) reads 

^ c[^j{pi +P3) • {P 2 +Pi)+c[f{pi +P3) ■P 2 +^ifG23{s,t,u)] 


with 


(16) 

(17) 


, (18) 


G23{s,t,u) = -92P2-PiiPl+P3)-{P2+P4) 

+2^3 [{Pl • P 2 ){Pl • P4)Pl ■ {P 2 + P4) + {pi Pa)] ■ (19) 

1“^ fsi 

The corresponding coefficients can be found in Table [2l The term survives only for 
inelastic scattering processes. 


3.2 One-loop contribution 

The one-loop connected graphs for D-cj) scattering are shown in the second line of Fig. 13.11 
All the vertices in the loop graphs originate from the Lagrangians T ^0 and which are 
free of unknown LECs. Similar to the tree-level amplitudes, it suffices to calculate the loop 
amplitudes for the 10 physical processes. All these loop amplitudes are listed in Appendix [B] 
which are expressed in terms of a set of one-loop integrals given in Appendix lAl 
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Table 1: The coefficients in the LO and NLO tree-level amplitudes of the 10 relevant physical 
processes. The Gell-Mann-Okubo mass relation, 3M^ = 4M|. — M^, is used to simplify the 
coefficients when necessary. 
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3.3 Renormalization 

In the previous sections, the IPI Feynman graphs are all calculated, which are related to 
the so-called amputated amplitudes. To derive the S'-matrix elements, one should perform 
wave function renormalization. Moreover, in the end, all the bare parameters should be 
replaced by the corresponding physical ones. 

3.3.1 Wave function renormalization 

To perform the wave function renormalization, one multiplies the external lines with the 
square roots of the wave function renormalization constants of the corresponding fields and 
takes them on the mass shell. In perturbation theory, if the calculation is done up to a certain 
order (up to 0{p^) in our case), the wave function renormalization is equivalent to taking 
the graphs in Fig. [2] into account. All the higher order contributions beyond the required 
accuracy are ignored. 



Figure 2: Feynman diagrams for the wave function renormalization at O {p^)- 
Hence, when taking wave function renormalization into consideration, the scattering am- 













Table 2: The coefficients in the NNLO tree-level amplitudes of the 10 relevant physical 
processes. 


physical process 
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plitude becomes 

A{s, t) = t) + t) + *) + Aolpi^^ + “^wf ' (^0) 

The first three terms are tree contribution given in Section 13.11 while the fourth term is the 
loop contribution discussed in Section [3^ and Appendix [Bj The last term A.wf('ST) corre¬ 
sponds to the contribution from the wave function renormalization. It can be obtained from 
the LO amplitude in combination with the wave function renormalization constants. For 
instance, considering the scattering process —)■ D 24 > 2 , it is given by 

= \{5Zd,+5Z^,+5Zd,+5Z^,)A^^^^{sA) , (21) 

with 5Z = Z — 1 and Z being the wave function renormalization constant up to the order 
considered. To be explicit, the wave function renormalization constants for D and Dg are 
Zd = Zds = 1 and for the Goldstone bosons are 


Z-n 


Zk 


Zn 


1 

1 

1 


7?2 

-^0 

1 

1 

•^0 


8L4[2M‘^ M^) + 8L^M^ '^K + > 

8Li^{2M‘^ + M^) + 8L^M]^ + 2 ^^ 4 ^*^ ’ 

8U{2Ml + Ml) + ^L5(4Mi - Ml) + Xk , 


( 22 ) 


where the tadpole loop integral Xi{i = vr, K, -q) can be found in Appendix [Al Note that in the 
above expressions, the ultraviolet (UV) divergence of the loop functions is not subtracted on 
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purpose. This is due to the fact that the 2^’s are not physical observables such that they might 
be divergent, namely the LECs L 4 and L 5 are not sufficient to absorb the UV divergence in 
those expressions. The UV divergence cancellation as well as the PCB terms absorption will 
be discussed in the following section at the level of the S'-matrix elements. As one will see, 
the S'-matrix elements are free of any divergence. 


3.3.2 Extended-On-Mass-Shell snbtraction scheme 


The loop integrals in the amplitude shown in Eq. (1201) is UV divergent, and we need 
renormalization to absorb the divergences by counterterms. Moreover, PCB terms show up in 
the chiral expansion if we use dimensional regularization with the MS scheme. It is necessary 
to get rid of them to have a good power counting. We will use the EOMS subtraction scheme 
which has the proper analyticity and correct power counting for the amplitudes. The essence 
of the EOMS scheme is to perform two subsequent subtractions: the MS subtraction and the 
EOMS finite subtraction. 

In the MS subtraction, the UV divergent parts are extracted and then cancelled by the 
divergences in the bare LECs, which are separated into finite and divergent parts as follows: 


where R = + 7 ^; — 1 — ln( 47 r) with 7 ^ being the Euler constant, and d is the space-time 

dimension. The coefficients at (f = 0, • • • , 5), f5j (j = 1,2, 3) and Pfc (A: = 4,5) are given by 


ao = 0, 


ai = 0 , 

/3i = 0, 


02 = 
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16’ 


/^s — 0, r4 
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Although the UV divergences have been removed so far, it is still not sufficient to get an 
amplitude that respects the power counting rule given by Eq. ([3]). The charmed mesons show 
up in the Feynman diagrams, say Loop(l) and Loop(2) in Fig. I3.lt generate the so-called 
PCB terms that often spoil the convergence of the chiral expansion [33]. The EOMS finite 
subtraction is used to get rid of those PCB terms. For each physical process given above, the 
PCB terms are easily obtained by replacing the loop function see Eq. (IB.ip . by 

namely Eq. (jC.3p . in the amplitudes and then performing the chiral expansion 
with respect to the small quantities. Note that the infrared regular parts for the required 
scalar loop integrals are also listed in Appendix O for easy reference. Eventually, the PCB 
terms are absorbed by decomposing the MS-renormalized LECs in the OijP') Lagrangian via 


id) — hi + 


167r2F2^° 


with the coefficient 5i (f = 0, • • • ,5) defined by 


-5o 


54 


5i = 0 , 52 = - 7 ^ + log ^ , 53 = 7 ^ - Tw log ^ , 
72 48 24 16 


35 7 , m2 ^ 35 7 , m2 

72M^ ^ 24M2 ® “ 48M2 “ 16M2 'JiF 


(24) 


(25) 


10 














The other LECs such as g^{n) and T^(/i) are untouched when performing the finite EOMS 
subtraction. 

After the two steps described above, we have obtained the full renormalized amplitudes. 
Eor the sake of easy practical usage, the chiral-limit D meson mass Mq and the chiral-limit 
decay constant Fq should be further related to the corresponding physical quantities according 
to the following expressions: 


M|, = Mo^ + 2 (ho + hi)M^ + 4 ho M| 


Mk 


— + 2 (ho — hi)M^ + 4 (ho -|- hi)M]^ 


= 




+ 1 : 


K 




Fn 


Fn 


(26) 

(27) 

(28) 


Here, we rewrite Fq in terms of rather than Fk and E^. This is the convention to be used 
throughout. Alternatively, one can also rewrites it in terms of F^ or F^^, and the difference 
is of higher order. The loop functions and LECs with a superscript r stand for their finite 
parts, namely, the contributions proportional to the UV divergence R are removed. 


4 Partial wave amplitudes and unitarization 

In this section, we will illustrate how to obtain partial wave amplitudes with definite 
strangeness S and isospin I from the 10 physical process amplitudes exhibited in the previous 
section in detail. Then we will discuss the unitarization of the scattering amplitudes using 
two different approaches. Based on the content of this section, it is straightforward to derive 
the S-wave scattering lengths, which will be discussed and compared with lattice data in the 
next section. 

4.1 Amplitudes for given strangeness and isospin 

The scattering amplitudes in the isospin basis can be classified by two quantum numbers, 
which are the strangeness S and isospin I of the scattering system. Hereafter, the scattering 
amplitudes with definite strangeness and isospin are called strangeness-isospin amplitudes 
for short. All the strangeness-isospin amplitudes can be related to the 10 amplitudes of the 
physical processes using crossing symmetry and isospin symmetry. 

We begin with the single-channel interactions. There are 4 single channels in total. The 
corresponding quantum numbers of {S,I) are (—1,0), (—1,1), (0,3/2) and (2,1/2). Their 
strangeness-isospin amplitudes are related to the physical-process amplitudes by 

= 2A£)+i^+^£)+i^+(tt,t,s) — n) , (29) 

= •^D0K-^D0K-is,t,u) , (30) 

= ■^D+^+^D+n+{s,t,u) , (31) 

Eor the coupled channels with {S,I) = (1,0), the strangeness-isospin amplitudes read 

= 2Aijox-^D°K-{'^^t^s) — Aij+x+^D+K+is,t,u) , (33) 
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For the coupled channels with {S,I) = (1,1), one has 


= Aj^+^o^^+^o{s,t,u) , (36) 

t, u) = •A£)+]^+^£)+j^+{s,t,u) , (37) 

•^D’K-^Ds7r(^’^’^) = \/2Ajj+j^-_^jyo„o(u,t,s) . (38) 

For (S,I) = (0,1/2), there are three channels: DTT,Dr] and DgK. The isospin relations are 
given by 


^(0,1/2) 

^( 0 , 1 / 2 ) 

'^Dri^Drj 

^( 0 , 1 / 2 ) 

•'^DsK^DsK 

4 ( 0 , 1 / 2 ) 

'^Drj^Dn 

4 ( 0 , 1 / 2 ) 

•'^DsK^Dit 

4 ( 0 , 1 / 2 ) 

■^DsK^Dri 


{s,t,u) 

is,t,u) 

is,t,u) 

(s,t,u) 

{s,t,u) 

{s,t,u) 


^Ad+tt+—>-D+tt+('^Aj^) 2'^D+tt+—>-D+tt+(^A:'^) , 

•^D+K+^D+K+('^AiS) , 

's/SAdo^^DO^o (s A:'^) , 


(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


4.2 Partial wave projection 

(SI) 

Each of the strangeness-isospin amplitudes can be denoted by {s,t)- Its partial 

wave projection with definite angular momentum i is given by 

Af’^\s)Di 4 >r^D 24>2 = ^ j cos 9 Pi(cos 9) jj^^^{s,t(s, cos 9)) . (45) 

Here, the Mandelstam variable t is expressed in terms of s and the scattering angle 9, 

t{s, cos 9) = ~ ^ ('^ ^D 2 ~ ^11) 

"a ’ "d,, Mfe) ■ (46) 

with A(a, b, c) = -|- 6^ -|- — 2ab — 26c — 2ac the Kallen function. From Eq. (j46h . one sees 

that at each of the thresholds of and he. when s takes one of the following two 

values 


Si = (Md^ -I- , S2 = (Md^ + M^j)^ , (47) 

t is independent of cos0. Taking s = si for instance, the 5-wave amplitude becomes 

This means that the 5-wave amplitude at threshold can be obtained directly from the full 
amplitude by setting the energy squared at its threshold value. However, note that this simple 
recipe can only be used for the single channel case. For coupled channels, it is necessary to 
perform the partial wave projection using Eq. (j45p . 
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Before ending this section, we remark that it is helpful to use matrix notation to denote the 
partial wave amplitudes with definite strangeness S and isospin I. In the matrix notation, the 
subscript —)• D24>2 is redundant. For single channels, this is apparent since the process 

is specified uniquely by {S,I). For coupled channels, taking {S,I) = (1,1) for example, there 
are four processes: DgTr —)• DgTT, DK —)• DK^ DK —)• Dgir and its time reversal process. 
Using time reversal invariance, one can write 

/ ^sjOsTT—^DsTT 
(s'joK^DsTr 

(SI) 

Later on, we will refer to the amplitudes for a given process in the isospin basis by A\ ’ {s)ij, 
with i and j being channel indices. Unitarization of the scattering amplitudes will be discussed 
in the matrix notation in the following. 




A^I’^\s)dK^Ds1t\ 

A^I’^\s)dK^DK j 


4.3 Unitarization 

Unitarization is often adopted to extend ChPT to higher energies. The unitarized am¬ 
plitudes sum up a series of s-channel loops, which correspond to the right-hand cut, and 
thus one would naively expect that they can be used for higher momenta as well as larger 
pion masses. Phenomenologically, it is now well-known that the unitarized amplitudes can 
well describe the scattering data for the pion and kaon systems up to 1.2 GeV, see, e.g.. 
Refs. [28l|44]. We thus expect that these amplitudes allows for a description of the lattice 
data at pion masses higher than the conventional ChPT. Yet, there is no rigorous proof a 
priori. For varying the quark masses (or equivalently the masses of the Goldstone bosons), it 
provides a way to performing the chiral extrapolation of lattice simulation results or studying 
the quark mass dependence of physical quantities. In the present work, we will consider two 
different versions of unitarization for the sake of comparison and for quantifying the inherent 
model-dependence of such approaches. For the sake of simplicity and generality, all the quan¬ 
tum number indices of the amplitudes such as S, I and I will be suppressed in this section. 
That is to say T, A, T and T, which will appear later on, are and 

respectively, for our case. 

The first approach we will use is the one proposed in Ref. [29], which is denoted by UChPT 
throughout this paper. In matrix form, the unitarized amplitude is given by 

r(s) = {i-r(s)-5(s)}-'-r(s), (50) 


where g{s) is a diagonal matrix g{s) = diag{g(s)i}, with i the channel index. The fundamental 
loop integral g{s)i reads 


g{s)i = i 


d^q _1_ 

(27r)4 (g2 _ Ml, + ie){{P - qf - M|, ie) ’ 


= 


(51) 


®Since the unitarization procedure is normally equivalent to a resummation of the scattering amplitudes in 
the s-channel, it breaks the crossing symmetry. Crossing symmetry can be restored using Roy-type equations, 
for an early attempt, see Ref. |43| . 
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Note that g{s)i is counted as 0{p) and its explicit expression is 


9{s)i 


1 

Idvr^ 


a(/r) + In 



• - Ml + Ml Ml 

2s Ml 


+ M|,, + Gi) — ln(—s + M|. — Mj). + cTj) 
+ ln(s + — Mq. + (Ti) — ln(—s — + M^. + at)] | , 


(52) 


with Gi = {[s — (M^. + — (M<^. — and p the renormalization scale. One 

can define a /r-independent parameter d = a(//) + ln(M|,y;U^), since a change of p in the 
logarithm can be compensated by a(/r). Notice that the parameter d in g(s) of Eq. (ISOl) 
cannot be absorbed by redefining the LECs. It is introduced through the dispersion integral 
along the right-hand cut, and is a free parameter in principle. The only constraint here is 
from the requirement of a proper power counting: while all other terms in Eq. (j52|) are of 
order O (p), d should be much smaller than 1 so that its presence will not cause a breaking of 
the power counting if we expand the resummed amplitude to a certain order, i.e. d = O (p). 
The kernel matrix T(s) can be obtained perturbatively by matching to the ChPT amplitudes 
order by order. Up to NNLO, it can be expressed as 

T(s) = -I- -I- • g{s) ■ (53) 

where ^^"^(s) (n = 1,2,3) stand for the partial wave amplitudes from the perturbative 
calculation with the superscript n denoting the chiral dimension. Notice that the right hand 
cut from the NNLO amplitude is subtracted in the last term in order to avoid double counting 
in the unitarization. In the function g{s) in the above equation, the subtraction constant d 
may be removed as it can be absorbed into the redefinition of the LECs in ^^^^(s). 

The other approach is the so-called inverse amplitude method (lAM) [401I411B1] . In our 
case, the lAM unitized amplitudes has the matrix form 


T{s) = f^^\s) • 


r(i)(s) -f(2)(s) 


fW(s) , 


(54) 


where 


f (1) (s) = (s) , (s) = (s) + (s) . (55) 

The above assignments guarantee that the unitarized amplitudes exactly obey unitarity when 
the perturbatively unitary equations are employed, i.e., 

ImA^^^(s)=0, ImA^^^(s)=0, Im A^^^(s) = A^^^(s)/i(s) A^^^(s)^ , (56) 

with p{s) = diag{p(s)i}, p{s)i = —qi/{8'Ks/s) and qt is the magnitude of the center-of-mass 
(CM) three-momentum in the channel. 


14 








5 Calculation of the scattering lengths 


5.1 Definition and pion mass dependence 

Given definite strangeness S and isospin I, the S-wave scattering lengths of the channel 
are related to the diagonal elements of the T-matrix, 0 

^ "87r(Mz), 

Here, Mjj. and denote the masses of the charmed meson and Goldstone boson cj) in 
the channel i, respectively, and (sth) stands for the 5-wave unitarized amplitude at 

threshold using either UGhPT given by Eq. (j5Up or lAM given by Eq. (I54p . 

Due to the short lifetime of the charmed meson, there are no experimental data for D-4> 
scattering lengths. Nevertheless, lattice QGD calculations in the last a few years provide very 
valuable information on the interaction between the charmed mesons and light pseudoscalar 
mesons [IMS]. Since the lattice calculations were performed at several unphysical pion 
masses, in order to describe these lattice data, one should know the pion mass dependence of 
the scattering lengths. This is achieved by replacing all the quantities in the expressions by 
the pion mass dependent ones. For the involved meson masses, we have 


n. - o o 

Mk = JmI + M2/2 , Md = Md + {hi + 2ho)^ , Mn^ = Md^ + 2ho^^ . (58) 

'' Md Md, 

Note that all the formulae shown above are of NLO for the pion mass dependence. Here the 
LEG hi can be fixed by the mass difference between D and Dg- Using these two equations, 
one has [12] 


hi = 


A{Ml - M2) 


= 0.4266 , 


(59) 


where the physical values for the meson masses are used, i.e., = 138 MeV, Mk = 496 MeV, 

Mo = 1867 MeV and Mo^ = 1968 MeVHxhe pion decay constant should also be substituted 
by [32] 

f 4/l/f2 8 472 1 

= Fo |l - 2^, - ^ [Ll{^i) + Ll{^i)] + , (60) 


where Mk is understood as the one in Eq. (1581) . and is a scale dependent function for the 
Goldstone boson (j) 




(61) 


^We are using the sign convention such that the scattering length for a repulsive interaction is negative. 

® In Eq. (I58II . although the formula we used for the kaon mass is a LO expression in SU(3) ChPT, it contains 
two parts: the part proportional to Bonis remains in the SU(2) chiral limit and is regarded as a LO 

contribution of the pion mass dependence, while the part related to Bom^/d ~ M^/2 vanishes in the SU(2) 
chiral limit and is thus a NLO contribution. In this sense, we spelled out the pion mass dependence for all of 
the masses and decay constants consistently up to the order M'^. 

®The mass of the r; is always expressed in terms of Mtt and Mk through the Gell-Mann-Okubo mass 
relation, 3M^ = 4M^ — M'^. 
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Figure 3: Chiral extrapolation of masses and decay constants. All the lattice data are obtained 
from the same ensembles, namely M007-M030. Data for Mk, and is taken from 
Ref. [IB] and the one for FV and Fk from Ref. [45]. Except for the data errors are so 
tiny that we do not show them explicitly in the plots. The vertical dashed line corresponds 
to the physical pion mass. 

So far, except for hi, the LECs L^, Lg and ho and the chiral limit quantities Mr, 
and Fq are all unknown. Since we will fit to the lattice results on the scattering lengths 
calculated in Ref. |16] . we choose to hx the above mentioned quantities from fitting to the 
lattice data calculated using the same gauge configurations. In addition, the kaon decay 
constant Fr data are also included to have a bigger data set for fixing Fq and 5 . The pion 
mass dependence of Fr is given by [32] 

Fr = Fo|i — ^ (Ftt + + fj,n) + ^2 + ^2 |) (62) 

with = (4M|. — Ml)I?) and Mr given by Eq. (l58l) . 

The lattice data for Mr, Md and are taken from Ref. [I 6 ]. There are four data sets 
for each quantity, corresponding to the four ensembles (labelled by M007, MOlO, M020 and 
M030) with pion mass approximately 301.1 MeV, 363.8 MeV, 511.0 MeV and 617.0 MeV, in 
order. Since the same ensembles are employed in Ref. [l5|, we take the data for and Jr 
from Ref. [45], where /tt = V^lBr and fR = V^Fr. Those lattice data are well described as 
shown in Eig. [3|0 when the parameters take the values given in Table [3| Our fitting values 
for L 45 are consistent with the determinations given in Refs. [46[l47j. Therein, the values are 
obtained at /r = M^, and the corresponding values transformed to /r = Mp can be found in 
Ref. dH. 

^°We have neglected the subtleties due to the use of mixed action gauge configurations in the lattice cal¬ 
culations, which in principle requires to use the partially quenched ChPT instead of the standard one for the 
chiral extrapolation, and the effect of finite lattice spacing, see Ref. [3^ . 
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Table 3: Parameters for chiral extrapolation. L 4 and are obtained at /z = Mp (= 
775.5 MeV). The masses and decay constant in the chiral limit are in units of MeV. /iq 
and hi and dimensionless. The asterisk marks an input value. 


Mk 

Md 

Md. 

ho 

hi 

Fo 

103 • LI 

103 • LI 

560.41 

1940.4 

2061.2 

0.0172 

0.4266* 

73.31 

0.0095 

1.3264 


5.2 Fits to lattice data on the scattering lengths 
5.2.1 Introduction to the fitting procedure 

Since all the necessary preparations are completed, we proceed to the description of the 
lattice QCD data of the S'-wave scattering lengths. There are two points to be discussed 
before carrying out the fits. 

The hrst one is related to the lattice data. From Ref. m, 20 data for 5 channels are 
available. Amongst the five channels, the DsTt with {S,I) = (1,1) can actually be coupled to 
the isovector DK channel while the other four are single channels. Although in Ref. m only 
the DgTT interpolating operator was constructed and used, the propagation of all the quarks 
should know about the presence of the coupled DK channel with (S',/) = (1,1) because the 
channel-coupling in this case does not require disconnected Wick contractions which were not 
included in Ref. [16j . Thus, we will describe the Dgir data using a coupled-channel unitarized 
amplitude. 

In addition, lattice QCD results were published in the last two years for two more channels: 
Dtt with (S,/) = (0,1/2) [IT] and DK with (S, /) = (1,0) [18]. These channels are more 
difficult since both of them involve disconnected Wick contractions, 0 but they are also more 
interesting as they can provide valuable information for the lightest scalar charmed mesons in 
the corresponding channels. The calculation for the Dtt scattering was performed using Nf = 
2 gauge configurations, and the DK calculation has results from both Nf = 2 and Nf = 2-\-l 
gauge configurations. Because the amplitudes derived here are based on SU(3) ChPT, we will 
only include in the fits the new result with Nf = 2-|-l, i.e. o!'^'k^dk ~ —1-33(20) fm obtained 
at Mtt = 156 MeV, and the unitarized amplitude used in the hts is obtained including the 
Dgi] coupled channel. Notice that these new lattice calculations use gauge configurations and 
actions different from those in Ref. m, the chiral-limit masses for the kaon and charmed 
mesons should take different values from those given in Table [3l Because the physical masses 
of the involved ground state mesons such as the kaon and charmed mesons were reproduced 
rather well with the lattice setup used in Ref. [I8| (for details, see Ref. |19]i. the chiral-limit 
values of the involved meson masses and Fq are determined by requiring them to coincide 
with the corresponding phyiscal values at the physical pion mass, namely, Mk = 486.3 MeV, 
Mf) = 1862.3 MeV, Mds = 1967.7 MeV and Fq = 76.23 MeV. The values for the LECs in 
the extrapolating expressions of these quantities are the same as those listed in Table O 

The other point concerns the LECs to be determined. There are 7 unknown LECs in 

is shown in Ref. [4^ that as long as the singly disconnected Wick contractions contribute, which is the 
case for the isoscalar DK channel, they are of LO in both the 1/A7 and chiral expansion. Therefore, they 
cannot be neglected. 
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total: /i 2 , / 13 , hi, / 15 , gi, §2 and 53 . As mentioned in the previous work m. /i 2 (/ 13 ) and 
hi (/ 15 ) are largely correlated. Therefore, redefinitions of the LECs are employed to reduce 
these correlations, which are 

/i24 = ^-2 + /14 , /135 = /13 + 2 /15 , hi = hiM'^ , /ig = h^M'^ . (63) 

The new parameters /i 24 , /i 35 , /i 4 and /ig will be determined in our fits. The average of the 
physical masses of the charmed D and Dg mesons, Md = is introduced 

to make the four new parameters dimensionless. Similarly, for the LECs from the NNLO 
contact terms, 52 and are largely correlated with each other, and it is better to redefine 
these LECs as 

923 = 92 - 2ff3 ! 9i = 9 iMd , 92 = 92Md , g's = gsM^ ■ (64) 

The parameters g[, 523 and g'^ have a dimension of inverse mass and will be fixed from 
fitting to the lattice data. One can fix g[ and 523 separately only when the coupled-channel 
unitarized amplitudes are used, i.e. from fitting to the lattice results of the Dgn and the 
isoscalar DK scattering lengths. The single-channel unitarized amplitudes is only sensitive 
to the combination 3123 = 523 — 9i, instead of g[ and g '23 separately, and g'^. 

5.2.2 Results 

We will try different fit procedures. In the fit UChPT- 6 (a), all of the 20 data points 
for 5 channels, with pion masses from 301 MeV up to 617 MeV, in Ref. m as well as the 
Nf = 2 -|- 1 datum for the isoscalar DK channel, with an almost physical pion mass of 
156 MeV, in Ref. [18] are taken into consideration. We notice that there are two possibilities 
for a scattering length to be negative in our sign convention: a repulsive interaction, and 
an attractive interaction with a bound state pole below the threshold. In the {S,I) = (1,0) 
channel, there is the well-known state i4*Q(2317) below the DK threshold which was not 
included as an explicit degree of freedom in our theory. Because the number of data is small 
but the number of parameters is large, a direct fit to these lattice data might result in solutions 
which are not physically acceptable. Eor instance, within the range of the parameters of a 
direct fit, the {S,I) = (1,0) DK channel could even be repulsive which is reflected by the 
fact that the kernel of the unitarized amplitude takes a positive value at the threshold. Given 
that the LO interaction in the corresponding DK channel is the most attractive one among 
all the charmed meson-Goldstone boson scattering processes, see Table II in Ref. [I 6 | for 
instance, we regard such a situation as unacceptable. Therefore, we put a constraint by hand 
requiring that when all the particles take their physical values there is a bound state pole 
in the {S,I) = (1,0) channel at 2317 MeV. Eollowing Ref. [T 6 |, this is done by adjusting 
the subtraction constant a in the loop function g{s) in the unitarized amplitude, Eq. (f50]i . 
to produce the pole at the right position. The resulting values of the LECs from the fit are 
shown in Table [H 

However, a pion mass larger than 600 MeV is definitely too large for the chiral extrapo¬ 
lation using the standard ChPT. The unitarized approach arguably has a larger convergence 
range than the standard ChPT. But the range is not known a priori. Therefore, for the sake 
of comparison, we perform another fit, denoted as UChPT- 6 (b), using the same method but 
excluding the lattice data at = 617 MeV. The fit results are shown in the third column 
of Table 01 One can see that the values of all the LECs from these two fits are similar, but 
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Figure 4: Comparison of the results of the 6-channel fits (without a prior x^) to the lattice 
data of the scattering lengths. UxPT-6(a): solid bine line with red band, UxPT-6(b): dashed 
blue line with green band. The filled circles are lattice results in Ref. m, and the filled sqnare 
(not included in the fits) and diamond are taken from Ref. |18] . 
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Figure 5: Comparison of the results of the 6-channel fits (with a prior to the lattice data 
of the scattering lengths. UxPT-6(a); solid blue line with red band, UxPT-6(b): dashed blue 
line with green band. The filled circles are lattice results in Ref. m, and the filled square 
(not included in the fits) and diamond are taken from Ref. [18] . 
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Table 4: Values of the LECs from the 6-channel fits using the method of UChPT. The /ij’s 
are dimensionless, and the 523 and g'3 are in GeV“^. 



UChPT-6(a) 

UChPT-6(b) 

UChPT-6(o') 

UChPT-6(6') 


no prior 

no prior 

with prior 

with prior 

^24 

n 70+0.10 
O’- ‘^-0.09 

f) 7f5+0-io 

O’- ‘O_0.09 

0-8^-ofo 

0.80+i° 

^35 

n 'yqH-0.50 
'^-0.38 

n oi +0.95 

0.43l°;i 

0.40l°;i 

K 

— 1 40+0-55 
-‘■•^^-0.57 

-1-561°:^^ 

1 qq+0-60 
-'-•^'^-0.60 

1 ^70+0.64 

' ^-0.63 

K 

— 11 47+2-24 

-15-38+®J 

—4 25'*'0’-05 

^■'0'4-0.87 

a'l 


—2 44+0’-57 

— 1 10+0-18 

-‘-•^'“’-0.35 

923 


^•O’O'-O.Sl 

—0 70+°-^® 

U./U_Q 24 

-1.481“;“ 

93 

9 1 9+0*55 
^•-‘-^-0.45 

9 qc: + 1-41 
^■°0>o.96 

0'-0’O-0.14 

0-58+?° 

xVd-o.f. 

31.52 _ 9 oc: 
21-7 ~ 

13.43 — 1 40 
16-7 ~ 

77.72-23.34 o oo 

21-7 ~ 0.00 

51.71-16.60 _ 3 90 


those from UChPT-6(b) have larger uncertainties as a result of being less constrained. The fit 
results from both hts are plotted in Fig.Sl The bands represent the variation of the scattering 
lengths with respect to the LECs within l-cr standard deviation. As we can see, both fits 
describe the lattice data reasonably well with the exception that the isoscalar DK scattering 
length around M.^ = 156 MeV is too large in comparison with the lattice result. However, 
both hts are consistent with the Nf = 2 lattice result for DK at a pion mass around 266 MeV 
which was not included in the hts. We notice that the lattice ensemble for the M-,^ = 156 MeV 
datum has a rather small volume with M^L ~ 2.3. It is a bit too small for Luscher’s hnite 
volume formalism to be strictly applicable, and thus this datum might bear a large system¬ 
atic uncertainty. The isospin-3/2 Dtt —)• Dtt scattering length vanishes at the chiral limit as 
required by chiral symmetry. Lattice discretization often breaks chiral symmetry. However, 
due to the use of the domain-wall action for the valence quarks in the lattice calculation of 
the pionic channels, the chiral behavior is protected in our case. For related discussions in 
mixed-action ChPT, we refer to Refs. |50H53j . 

In both fits, the values of all the LECs except for /15 turn out to be of a natural size. 
However, the absolute value of the dimensionless LEG h'^ is too large to be natural. This 
means that the absolute value of h'^ is so large that this single term would give a contribution 
larger than the LO amplitude. It would spoil the convergence, and thus the perturbative 
expansion, at least for some quantities (although for some other quantities, due to hne-tuned 
cancellation the sum of the NLO contribution could still be much smaller than the LO one). 
Therefore, we try to constrain all the LECs to natural values following Ref. |54j which discusses 
the use of the Bayesian method in effective held theories. Following that paper, the so-called 
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augmented chi-squared can be defined by@ 

xLg = x’ + Xpriot , (65) 

where is the usual chi-squared used in the standard least chi-squared fit and Xprior i® ^ 
prior chi-squared encoding the naturalness requirement of the fit parameters. In our specific 
case, the XpOor i® ®®i' ®^™ squares of the fit LECs. This means that we require the 

dimensionless LECs hi '^’s to be 0(1) and g'^s to be 0(1 GeV~ ). The results by minimizing 
the augmented chi-squared are listed in the last two columns in Table 01 denoted as UChPT- 
6(a') and UChPT-6(6'), where the values for are given with Xprior subtracted. One sees 
that the value of /ig gets more natural at the price of a larger A comparison of the 
scattering lengths with the lattice data in various channels is given in Fig. 01 and one can see 
that the lattice data in all six channels can still be described reasonably well. 

It turns out that in all of these fits \h'^\ > which is consistent with the N^. counting 
1^4! = O (\h'^\/Nc) [H]. The values of the hj’s are different from those obtained in Ref. [T6] . 
The reason may be attributed to the use of the EOMS scheme in this work, and all of 112,3,4,5 
absorb a power counting breaking contribution, see Eq. (1241) . For the case of the UsTt, the 
scattering length does not vanish at the limit of a vanishing pion mass. This is due to 
the presence of the L^RT-loop in the coupled-channel amplitude which has a nonvanishing 
contribution in the SU(2) chiral limit. We have checked that the elastic contribution tends 
to zero as M-j^ approaches zero as required by chiral symmetry. 


Table 5: Values of the LECs from the 4-channel fits using both the methods of UChPT and 
lAM. The /ij’s are dimensionless, and the (7123 = 523 ~ g' and g'^ are in GeV“^. 


UChPT-4 IAM-4 


^24 

'^•^'-’-0.10 

Q co+0.07 

^35 


-0.591};°^ 


1 90 + 1.03 
-'-•^'^-1.08 

0.6+°:“ 

K 

-3.09+^;^^ 

-6.08l®;05 

9123 

n 10+0.18 

U.Z,0_q22 

CO 

-'-•'^-‘■-0.86 

1 42+1-0® 
-'■•+^-1.10 


xVd.o.f. tS = l-36 t^ = l-40 


For comparison, we also perform fits with just the four single-channel data, i.e. the 
Dgir and isoscalar DK data are excluded. For this case, we use two different unitarization 
methods: UChPT, to be denoted as UChPT-4, and lAM, to be denoted as IAM-4. We did not 
use the lAM approach in the 6-channel fits because this approach is not suitable to unitarize 

^^The method in Ref. was only derived for the case that the dependence on the parameters to be fitted 
is linear. Although our case is non-linear and thus the augmented lacks a strict statistical meaning, we still 
try this method as the defined in this way comprises a ’’naturalness prior” so as to favor natural values for 
the LECs. 
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Figure 6: Comparison of the results of the 4-channel fits to the lattice data of the scattering 
lengths. UxPT-4: solid red line with blue band, IAM-4: dashed red line with green band. 
The lattice data are taken from Ref. m¬ 


a perturbative amplitude with a zero LO contribution. As can be seen from Eq. ()54l) . if the 
LO amplitude vanishes the unitarized one will vanish as well. This happens to the case of 
the DgT^. The UChPT approach is free of this problem. The results of these two fits are 
compiled in Table |5l Notice that in this case g'^ and (723 cannot be determined separately, 
and the effective combined parameter is <7123 = 523 ~ g'l ■ One sees that the values of LECs 
from the hts using different unitarization methods are consistent with each other, El but are 
only marginally consistent with those in the 6-channel fits. In addition, the uncertainties are 
quite large. More lattice simulations are apparently necessary to pin down the LEG values. A 
comparison of the results of the 4-channel fits to the lattice data in these channels are plotted 
in Fig. m 

For reference, the values for the scattering lengths extrapolated to the physical pion mass 
are presented in Table [U The chiral limit values in Table [3| are adopted for all the 16 channels 
when performing the chiral extrapolation . Here we only show the results using the 6-channel 
fits to the data with the pion mass up to 511 MeV, i.e. UChPT-6(b) and UChPT-6(6'). We 
notice that the numerical results of the scattering lengths extrapolated to the physical pion 

^^However, not all of the LECs in these different unitarization methods ought to take the same values. One 
can see this by expanding the lAM resummed amplitude up to O (p®). Considering the single channel case 
for simplicity, one has TiamIs) = -I- -f -I- (s) + O (p*)- It is different from 

that of UChPT, TuchPTfs) = -I -+ 0 Thus, the LECs in the O (p®) Lagrangian 

could take different values. 
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Table 6: Predictions of the scattering lengths at physical pion mass using the LECs determined 
in the 6-channel fits UChPT-6(b) and UChPT-6(6') in units of fm. 



UChPT-6(b) 

UChPT-6(6') 

(-1.0) 

‘^DK^DK 

'11-0.31 

n 0-1-1-0.15 
n-ii’i-o.is 

(-1.1) 

‘^DK^DK 


11 -^11-0.02 

(0,1) 


n 42-I-0-04 

u-^^-o.os 

„(0.|) 

^Dri^Dri 

—0 18+0-0^ - 1 - 9 0 OO'*'^'*^^ 
'-'•i°-o.04 '■u-un-o.oo 

-0-2i;n!±f0.0ll0;0J 

(o.|) 

^DsK^DsK 

-‘■•'^'-0.04 ' '-'■111-0.02 

-0.471°;°? ±9 0.5010:18 

(0,f) 

—0 lA+O-lll 

H'l^-O.Ol 

-0 1 5+0-°l 
H-lH-O.Ol 

(1.0) 

^DK^DK 

1-11^-0.03 

l.OU-o 26 

(1.0) 

^Dstj^DsTI 

—0 ±90 ni+ 0 -'ii 

u-n^-o.os 'u-ui-o.oo 

—0 7R+'1-'13 _i_ n OK-I-O.oo 

H-'11-0.05 ± * U-UH-O.Ol 

(1.1) 

^DsTV^DsTT 

n-ni-o.oi 

—0 ni+o-oi 
U-Ul-O.Ol 

(1,1) 

‘^DK^DK 

— 1 ll+'l-23 _|_in ' 77 - 1 - 0.27 

—0 82'’'°'^° ± 9 1 64'''° °^ 
ll-i’^-0.38 ^ ' l-H^-O.ll 

'^DsK^DsK 

11 -^11-0.02 

—0 -^o+O-Oi 


masses in some channels differ from those obtained in Ref. m- This could indicate that 
the uncertainties are underestimated as the SU(3) formalism for UChPT was applied to pion 
masses higher than 500 MeV. We expect that the situation will improve when lattice results 
at lower pion masses are availble. 

5.2.3 Contribution of vector charmed mesons 

In this section, contributions from vector charmed mesons will be included explicitly in 
order to quantify their influences on the S-wave scattering lengths. The diagrams that survive 
in the heavy quark limit, see also Ref. [23], are taken into account and shown in Fig. 15.2.31 
Those diagrams vanishing in the heavy quark limit are suppressed by 1 /rric and therefore are 
neglected. We denote the vector charmed mesons by D* = (D*^, D*~^ , and the vertices 

involved in Fig. 15.2.31 are described by the following Lagrangian, 

Cd*d^ = + i g - D , ( 66 ) 

where the covariant derivatives acting on D* are analogous to those defined in Eq. m- 
Further, Mg is the mass of D* in the chiral limit. The relation between the axial coupling 
constant g defined here and the coupling g which is employed usually in the heavy meson 
ChPT [7I [9l[55] is g = y/M^Mj^g. Following Ref. |56|, we take g = 0.570±0.006, determined 
by calculating the decay width of the process D*~^ —)• D^7r~^, and then one gets g ~ (1103.3 ± 
11.6) MeV. The calculations of the Feynman diagrams in Fig. 15.2.31 are straightforward but 
the analytical results are too lengthy to be shown here. Similar to Eq. (|58p . the pion-mass 
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dependence of the D* and D* masses reads 


o _ o — jy± 

Md* = Mo* + {hi + 2ho) o ^ , Mog* = Mo^* + 2ho——— , (67) 

Mo* Mo,* 

where /iq and hi are the analogues of /iq and hi, respectively. In the heavy quark limit, one has 
hi = hi and ho = ho. As discussed in Ref. [26], the breaking of heavy quark spin symmetry is 
only about 3%. Therefore, to a good approximation, we impose these two heavy-quark limit 
relations. The masses of the vector charmed mesons in the limit of —)• 0, i.e. Mo* and 

Mo,*, are related to the corresponding ones of the pseudoscalar charmed mesons via 

Mo*^^^- - Mo^^^- ^ Mo* - Mo , Mo,*^"^^- - Mo,^^^' ^ Mo,* - Mo, , ( 68 ) 

with Mo*^^^' = 2008.6 MeV and Mo*^^^' = 2112.3 MeV, which are physical masses for D* 
and D*, respectively. In parallel to the four kinds of 6-channel fits in the previous section, 
we refit the 5-wave scattering lengths and the results are shown in Table [T] In each case, the 
LECs as well as the chi-squared are almost same as before. This implies that the influence 
of D* to the 5-wave scattering lengths is marginal and it is a good approximation to exclude 
them in the calculation. 

Table 7: Values of the LECs from the 6-channel fits (including explicit D*) using the method 
of UChPT. The hj’s are dimensionless, and the g'l, 523 and g'^ are in GeV“^. 



UChPT-6(a) 

UChPT-6(b) 

UChPT-6(a') 

UChPT-6(6') 


no prior 

no prior 

with prior 

with prior 

^24 

6 -o6_o.o8 

0.801°:^° 

0.851°;“ 

0.85l°;“ 

hsb 

n Q 9 + 0 . 6 O 

0.98l°|I 

0.50l°;i 

0.59l°;°° 

h'^ 

1 o'7H-0.52 
-0.51 

-‘-•^^-0.62 

1 99+0.58 
-'-•^^-0.57 

-1-59+^? 

K 

— 11 f) 1 + 2-53 

-15.06l^J? 

q oy+0.67 
-0.69 

'^•48_o.83 

a'l 

^■^^-0.36 

—2 

^•'^a_o 61 

— 1 411 + 0-20 

i-yo_o 43 

523 

— 2 

^•^'J-0.31 

Q qoH- 0.46 
^•^o_o_48 

— 1 10 '*'°'^^ 
■'-•■'-U-0.31 

1 c:i-H0.40 
-'■•^■'--0.45 

53 

9 1 c:-H0.60 
^•-‘■^-0.49 


0 qi+015 

'3-a-‘--o.i5 

U-aD o.l9 

xVd.o.f. 

29.36 — 2 10 
21-7 “ 

13.75 _ 1 crq 
16-7 ~ 

74.22-2L56 _ g 

50.06-15.96 0 

16-7 —0.0) 


6 Summary and outlook 

We have computed the D-cj) scattering amplitude that is valid up to the NNLO in the 
chiral expansion within the framework of ChPT. The complete analytical expressions for the 
amplitudes are given using a renormalization procedure with the EOMS subtraction scheme. 
We show explicitly that the UV divergences and the PCB terms, both of which stem from 
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( 6 ) 


(d) 

Figure 7: Feynman diagrams (including the vector charmed mesons) that survive in the heavy 
quark limit. 

the loops, can be absorbed into the LECs. We then obtained the EOMS-renormalized D-cj) 
scattering amplitudes which are independent of the renormalization scale and possess good 
properties such as correct power counting and proper analyticity. 

In order to describe the lattice data on the S-wave scattering lengths at relatively high pion 
masses and to account for the nonperturbative nature in the channels like the {S,I) = (1,0) 
DK, the aforementioned perturbative amplitudes are inserted into a unitarization procedure 
to perform the chiral extrapolation from large unphysical light quark masses down to the 
SU(2) chiral limit. We tried different fitting procedures with and without a naturalness con¬ 
straint. It turns out that the absolute value of h'^ could be quite large if the naturalness 
constraint is not put by hand. We want to stress that more lattice simulations in different 
channels are necessary for a better determination of the involved LECs and a better under¬ 
standing of the scalar and axial-vector charmed mesons. When the LECs are well constrained, 
we can make reliable predictions in the channels which have not been calculated on the lattice 
and in the bottom sector utilizing heavy quark spin symmetry. 



Acknowledgements 

We would like to thank C. Hanhart and Tom Luu for useful discussions and comments. 
This work is supported in part by DFG and NSFC through funds provided to the Sino- 
German CRC 110 “Symmetries and the Emergence of Structure in QCD” (NSFC Grant No. 
11261130311), and by NSFG (Grant No. 11165005). The work of UGM was supported in part 
by The Chinese Academy of Sciences (CAS) President’s International Fellowship Initiative 
(PIFI) grant no. 2015VMA076. 


A Definition of one-loop integrals 

In this appendix, all the relevant one-loop integrals are defined. For the current case, 
only one- and two-point loop functions are involved. As is well known, each tensor one-loop 
integral can be expressed as a linear sum of scalar one-loop integrals by using the method of 
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Passarino-Veltmann (PV) decomposition [57]. Hence, if the explicit expressions of the scalar 
one-loop integrals are known, the loop amplitudes can be obtained analytically. 

Throughout this work, the ultraviolet divergence is contained in the quantity R which is 
defined by 

+ 7 e - 1 - ln(47r) , (A.l) 

with 7 ^ the Euler constant and d the space-time dimension. In addition, we will denote the 
renormalization scale by fi. In terms of these notations, various loop integrals involved in the 
calculations are given as follows: 

• One-point loop function: 

r ^ _ 1 _ 

i J {2 'k)’^ k'^ — -|- i0+ IGvr^ 


R + \n^ 


(A.2) 


• Two-point loop function for unequal masses: (M^ > Mb) 

~ ~r J (27r)'^ (fcs - M2 + i0+) [{k + p)2 - m2 + z0+] ’ 

where the PV coefficients are given by 
Kbip^) = 

Ktip^) = 

= 

where we have defined Aab = and Tjab = M^ + . The scalar two-point 

one-loop function Rab{p^) has the following analytical form. 


^ \Xa -1b- {p^ + Aab)Rab{p‘^)] , 


1 


12p2 

1 1 
~ 167r2 18 


{ {p^ -h Aab)1a + (jp^ - Aab)1b + [dp^ - (j? + A^fe)^] Rab{p^)} 


{p^ - , 


{-(p2 + Aab)1a + (2/ + Aab)1b " [p^ - {p^ + Aab)^] Rabip^)} 


1 

+ . „ n . „ n “ ^1‘ab) ) 


167r2 18p2 


with 


T-iabiP^) 


1 


IGvr^ 



-i? + 1 - In ^ + 
P 


.Ml 

2p2 m2 



Mfe)2 


Pafe(p^)ln 


Pab{p^) - 1 ' 
PabijP^) + 1 _ 


Pab{P^) = 


p2 - [Mg + Mfc)2 

p^-iMa-Mb)^' 


(A.4) 


(A.5) 


To get the imaginary part above the threshold properly, one should take the branch cut 
for the logarithm along the negative real axis. 
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Two-point loop function for equal masses: 




4-d 


dU 




i J {27:)d {k^ - + i0+) [{k - pY - + i0+] ’ 

where the Passarino-Veltmann coefficients are given by 

“ P^)MP) + Ta + - P^) • 


1 1 


^ Rp" - +IJ + 15^2 I8p2 


(p" - 6MS , 


In this case, the scalar two-point one-loop function Ja{p^) a much simpler analytical 
form, 


JaiP^) = 


1 


IGvr^ 


-R+1-\ti^ + (Ta(p^) In I I 

O-a(P^) + IJ 


2 / 4M2 

> 0-a(p) = Wl- 




(A.6) 


B Loop amplitudes without explicit charmed vector mensons 

In order to express the loop amplitude in a short form, the following abbreviation is 
adopted, 

(s,i) = [3(s - Ml) + {s- Ml)] Xd-{s- T.,,)‘^ncd{s) + 2 (t - 2MI) U^is) 

+2 (s - (s - Sfce) Hi,{s) -{s- Aat)^ nllis) . (B.l) 

We first list the loop amplitudes concerning the elastic scattering processes. 

• D°K- D°K- 


+(s - u) {Xr, + 2 Xk + x^) - 4(s - u) [Jl\t) + 2J^\t)] |. 


(B.2) 


• D+K+ D+K+ 

- 4(s - u) [Jl\t) - J^\t)] } . 

(B.3) 

• T>+7r+ T)+7r+ 

(S> i) = *) 

+ ^(s - (2X, +X^) - 4(s - u) [2jl\t) + J^\t)] |.(B.4) 
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.loop 

'^D+ri^D+rj 


{s,t) 


1 

16^4 


3 (DsK) 

2 -^ Dr, 


(S) t) + -F ^ 


{D,K) 

Dr, 


{u,t) 


(B.5) 


D+K+ DtK' 




16F4 




+(s — ti) (Xjy + 2X/f +Xn-) — 12(s — u)J^^ it) 


(B.6) 


^ -0+ r] 




D+vrO ^ D+vrO 


4 loop 

D+ttO^.D+ttOI 


1 




(DK), 


(B.7) 


(B.8) 


As for the inelastic processes, the amplitudes become a little more complicated. To reduce 
them, we further need 

^abcdi^’'^) = 2 ^bd^ef{s) + - (Aac “ AM)^'Hg]-(s) - A^ (Aac “ ^bd)F-lf{s) . (B.9) 

In the above equation, the letters a and b (c and d) label the incoming (outgoing) particles, 
while e and / mark the particles in the loop. This convention also holds for the abbreviations 
^^ab,cd(^^ t) and ,t), whose explicit expressions are given by 


= ^ac <! + if [Hefit) + nlf{t)] - - tnl}{t) \, 


‘'ab^cd^ 

-(e/) r 
''ab,cd\ 


= -3(s - u)'Hgj(t) - Aq 


2 

,00, 


-{GJ^de + QAdf-l3Abd)nlf{t) 


+ -{3T,de + 3Adf — 2Abd)'Hef{t) — 3AbdHef{t) 


(B.IO) 


In combination with the abovementioned notations, the inelastic one-loop scattering ampli¬ 
tudes are given as follows: 

• D^rj —)• D^tt^ 


■^^doI^do-kO ('®’ ('®> ('5> t) +2 t)F{s^ u) 


(B.ll) 
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16F4 \ 2 
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-(‘cSSD.(».‘)+^':gSD.(».‘)) 


D+iv:- ^ 




(B.12) 


{1 ’{<-. *) + <■">(». ‘) + 2eSS,(». *)) 

+ - '^^DsK,Dr^isA)) 

\ '^DsK,dA°^^> ^ '^DsK,dA^^^'J 
^(5M2 + 8Mi - m 2) (2^],,(t) + 


+- 


A 


Ad, 


^(m 2 - 4Mi + m 2) (2^i^(t) + n,K{t)) 


+ 


(X, + X^ - 2Ik) - 4^ (X^ + x^ + 6Ik) 

o 4 


(B.13) 


C Infrared regular parts of the loop integrals 

The following expressions for the infrared regular parts are taken from Ref. |58] with the 
nucleon mass (pion) mass) replaced by the D meson (Goldstone boson) mass; 

• one-point: a € {D, Dg} 


^reg. 


AP Af2 

—2_in—2- 

Ibvr^ 


• two-point; a G {D, Dg} and b G {vr, K, r]} 


(C.l) 


Kf (-) = 


1 


IGvr^ 
1 


^ , M2\ s-M2 1 


2M2 167r2 


m2 

1-log^ 

fJ- 


+ 


327r2 




m2 


M2 


m2 

log^ 

7X2 


+ 0(p3). (C.2) 
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The power counting breaking term of is of 0{p^) and its explicit form reads 

1 




IGvr^ 
-{s-Mcf 


2{s - M^){s - M^) 


A'P 

1-log^ 


'll ; 

2{t- 2Mi)M^ 


-is- Mlf 


111 

"9 6 ^J? 


1 m2 
m2 


(C.3) 


Since the difference between M2 and is at least Oip^), the above expression can be 
reduced to a simpler form 




1 


1447r2 


2\2' 


2 (t - 2M^) m2 - 35 (s - m2) 


+3 


m2 


7 is-Ml) -it-2Ml)Ml log^ 
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(C.4) 
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